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Abstract: We study a phenomenological dark energy model which is rooted in the 
Veneziano ghost of QCD. In this dark energy model, the energy density of dark energy 
is proportional to Hubble parameter and the proportional coefficient is of the order Aq CD , 
where Aqcd is the mass scale of QCD. The universe has a de Sitter phase at late time 
and begins to accelerate at redshift around z acc ~ 0.6. We also fit this model and give the 
constraints on model parameters, with current observational data including Snla, BAO, 
CMB, BBN and Hubble parameter data. We find that the squared sound speed of the 
dark energy is negative, which may cause an instability. We also study the cosmological 
evolution of the dark energy with interaction with cold dark matter. 
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1. Introduction 

It has been more than ten years since our universe was found to be in accelerating expansion 
Q. A new energy component of the universe, called dark energy (DE), is needed to explain 
this acceleration. The simplest model of DE is the cosmological constant, which is a key 
ingredient in the ACDM model. Although the ACDM model is consistent very well with 
all observational data, it faces with the fine tuning problem Q|. Plenty of other DE models 
have also been proposed g, g, |, §, [7|, [5| g| , but almost all of them expl ain the acceleration 
expansion either by introducing new degree(s) of freedom or by modifying gravity. 

Recently a new DE model, so-called Veneziano ghost DE, has been proposed || ||. 
The key ingredient of this new model is that the Veneziano ghost, which is unphysical 
in the usual Minkowski spacetime QFT, exhibits important physical effects in dynamical 
spacetime or spacetime with non-trivial topology. Veneziano ghost is supposed to exist for 
solving the U(l) problem in low-energy effective theory of QCD |ll], [l3|, [L4|. The 



ghost has no contribution to the vacuum energy density in Minkowski spacetime, but in 
curved spacetime it gives rise to a small vacuum energy density proportional to Aq CD H, 
where Aqcd is QCD mass scale and H is Hubble parameter. This small vacuum energy 
density expects to play some role in the evolution of the universe. Because this model is 
totally embedded in standard model and general relativity, one needs not to introduce any 
new parameter, new degree of freedom or to modify gravity. With Aqcd ~ WOMeV and 
H ~ 10~ 33 ey, Aq CD H gives the right order of observed DE energy density. This numerical 
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coincidence is impressive and also means that this model gets rid of fine tuning problem 
B B|. Actually, several authors have already suggested DE model with energy density 
proportional to h? QCD H in different physical contexts |1| [U], ||, |0|, ||l|, ||, |0] 
such as QCD trace anomaly, gluon condensate of quantum chromodynamics, modifying 
gravity and so on. 

In this work, we investigate the phenomenological model with energy density of DE 
Pde proportional to Hubble parameter H. We study the cosmological evolution of the DE 
model with/without interaction between the DE and dark matter. We analytically and 
numerically compute some quantities such as scale factor a, pde, squared adiabatic speed 
of sound Cg and so on. Also we fit this model with current observational data and give 
constraints on the model parameters. 

The note is organized as follows. In section 2 we study the dynamical evolution of 
the DE model. In section 3, we fit this model with current observational data and discuss 
the fitting results. The data used are Union II Snla sample [22], BAO data from SDSS 
DR7 H|, CMB data (R, l a , z* 
and Big Bang Nucleosynthesis (BBN) (27 



from WMAP7 ||, 12 Hubble evolution data |6| 
In section 4, we calculate <? s and find it is 



always negative. Negative c 2 s may give rise to the instability problem. On the one hand, 
to try to avoid this problem, on the other hand, to further study the ghost DE model, 
we introduce the interaction between DE and cold dark matter (CDM), and study the 
dynamical evolution in this case. We summary our work and give some discussions in 
section 5. 



2. Dynamics of ghost dark energy 

To study the dynamics of the DE model, we consider a flat FRW universe with only two 
energy components, CDM and DE and neglect radiation and baryon temporarily in this 
section. We will include the radiation and baryon when fit the model with observational 
data in section 4. 

In this ghost DE model, the energy density of DE is given by pde = otH, where a is a 
constant with dimension [energy] 3 , and roughly of order of Aq^ where Kqcd ~ lOOMeV 
is QCD mass scale. Arming with this DE density, the Friedman equation reads 

H 2 = ^-(aH + p m ), (2.1) 

where p m is energy density of CDM, whose continuity equation gives 

p m + 3Hp m = => p m = pmoa~ 3 . (2.2) 

We have set ao = 1 and the subscript stands for the present value of some quantity. From 
(|2.1|) and (p.2[), we can obtain the Raychaudhuri equation 



H + H 
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(2.3) 



- 2 - 



Solve the Friedman equation, we have 




H± = -^—a ± a -—a + -— p m0 a- 3 - (2.4) 



There are two branches, H + represents an expansion solution, while iT_ a contraction one. 
We neglect the latter since it goes against the observation, and for simplicity, write H + as 
H in what follows. 

To facilitate our discussion, we define a characteristic scale factor a* 

1 n I 

3 
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MeV c 

2~ 



where we have taken Hq = 10 _33 ey, Q m o and Qdeo are the dimensionless energy density 
of CDM and DE, respectively. One can see shortly that actually a* is the transition point 
when the universe transits from the dust phase to a de Sitter phase. If we assume £~l m o = \ 
and Qdeo = j, by definition, we get roughly o* ~ 1, which means that the transition occurs 
just at present. Therefore, we will take a* = 1 throughout this section. And especially 
a ~ (lOMeV) 3 if a* ~ 1. 

At early epoch a < a» ~ 1 , the a~ 3 term dominates in ( |2.4[) , the Hubble parameter 

3 

behaves like H ~ which means that the universe is in a dust phase. While at late 
epoch a ^> a* ~ 1, the a term dominates in (12.4) , as a result, = const, which says that 
the universe enters a de Sitter phase at late time. Here a* is the transition point between 
these two phases as mentioned above. 



We can solve (2.4) analytically as 



3 



4irGa(t - U) = -x 3 + x 3 \A + x~ 3 + -lnx + ln (l + Vl + 



where x = ^- and is the initial time when a(ij) = 0. At early time id, AnGa (t — t») = 
2x2; while at late time i>1, AirGa (t — tj) = |lnx. These asymptotic behaviors agree 
with the previous argument. The numerical relation t ~ a is plotted in Figure [l]. 
Using Q2.4|) and the definition of puEi the energy density of DE is 



rr 47TG 2 

Pde = oil = — - — a 



, a- n 3 
1 + 1/1+ . 

a 



The behavior of /9d_e in terms of a is shown in Figure |2[ The equation of state (EoS) of 
the ghost DE is given by 

"Hi^- 1 -^' 1^=^ — ^ — ^l- 1 < 2 - 5 ' 

J —5 a < a* 
1—1 a » a* 
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Figure 1: A-KGa(t — U) ~ a, where a* = 1 




From the asymptotic behavior of w, we can see that the DE acts like a cosmological constant 
at late time. We plot the relation u) ~ a in Figure |3| From the figure, we can see that u 
can never cross —1, which is similar to the behavior of quintessence. We will show that this 
behavior can be altered in the presence of interaction between DE and CDM in section 4. 
ui varies from — ^ at early time to —1 at late time, which is similar to freezing quintessence 
model [p9| . We can also find that oo has a sharp variation round a = a*. It is easy to 
understand if we rewrite the expression of to as 

*+»> = -i^ = -# = (*-'>■■ <*» 

from this equation, we can see that in this model the EOS of DE tightly relates to the 
variation of Hubble parameter, which is quite different in different phases of the universe. 
For instance, in the dust phase, (-f^ -1 ) ~ f ; while in the de Sitter phase, (H~ l ) ~ 0. 
Therefore, there will be a jump from — ^ to —1 when the universe transits from the dust 
phase to the de Sitter phase. 

One can easily show that the EoS of the ghost DE model is u = — ^ . x . The value 
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Figure 3: u ~ a, where a* = 1 



wo at present is 

W ° (a = 1} = o 1 2 = ~n~~TT' (2 ' 7) 
where we have used a* = ( 4 i^ m " ) 3 in the first equality. This relation is important and 

V DEO / 

helpful to understand the fitting results in section 3. 

Also from Raychaudhuri equation ( [2.3D , we can get the total equation of state of the 
universe 

1 2 R 1 Ju 9 / a < a * (O Z\ 

u tot = -I- -—^ = 1 + 2^ = \ 2.8) 
3 H z 1—1 a > a*. 

wtot decreases monotonically from to —1, which means that the expansion of the universe 
switches from deceleration at early epoch to acceleration at late epoch. The conversion 
occurs at a acc when u tot (a acc ) = — ^. From (^) and fl2.8| ), we can calculate a acc 
which means that the universe begin to accelerate at redshift z acc ~ 1. 



3. Data Fitting 
3.1 Model 

In order to fit the model with current observational data, we consider a more realistic model 
which includes DE, CDM, radiation and baryon in a flat FRW universe in this section. In 
this case, the Friedman equation reads 

rrl 87TG , rr 

H = — — [aH + p DM + Pb + Pr) , 

which can be rewritten as 

H 1 Al 

e = — = -n DE0 + y -n 2 DE0 + (o DM0 + fiw) (i + zf + n r0 (i + z) 4 , 
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where Qdmo, ^&o> ^t-o are present values of dimensionless energy density for CDM, baryon 
and radiation, respectively. Qdeo = 8 3# a * s dimensionless energy density of DE at present. 
Energy density of baryon and CDM are always written together as £Idmo + ^60 = ^mo- 
Notice that Qdeo + ^ m o + ^ro = 1 since we assume a flat universe. The energy density of 
radiation is the sum of those of photons and relativistic neutrinos 

O r0 = ^ 7 o(l + 0.22717V n ), 

where N n = 3.04 is the effective number of neutrino species and f2 7 o = 2.469 x 10~ 5 /i~ 2 
for T cmb = 2.725K (h = H /100 Mpc ■ km ■ s' 1 ). 

It is worth noticing that from the definition of dimensionless energy density of DE and 



flatness of our universe we can get an important relation (see also (4.4)) 



8irGa 

(1 - il m0 ) H = — - — = const, (3.1) 

where we have neglected Q r o, which is very small compared to O m o. It means that pa- 
rameters O m o, h and a are closely related, a can be expressed in terms of f2 m o and h. We 
will choose £l m o, h and Qbo as free parameters of the model in the following analysis. This 
relation also infers that there exists a strong degeneracy between £l m Q and h, as shown in 
subsection 3.3. 

3.2 Sets of Observational data 

We fit our model by employing some observational data including Snla, BAO, CMB, Hubble 
evolution data and BBN. 

The data for Snla are the 557 Uion II sample [p2| . Xsn f° r Snla is obtained by 
comparing theoretical distance modulus [ith{z) = 51og 10 [(l + z) dx/E{x)] + fiQ (/xo = 
42.384 — 51og 10 h) with observed fi i, of supernovae: 

2 _ ^ \Pth(zi) - ^ob(z-i)} 2 



a 2 (zi 



To reduce the effect of fiQ, we expand Xsn with respect to 1 3C ] : 

X 2 sn = A + 25/io + Cn 2 (3.2) 



where 



y-v [Vthjzi ] = 0) - flob(Zi)} 2 



CT 2 (Z 1 ) 

IHh(Zi',t4) = 0) - Vob(Zi) 



CT 2 (Zi) 



|3.2|) has a minimum as 



~2 _ 2 __ a _ d2 ip 
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which is independent of /io- In fact, it is equivalent to performing an uniform marginaliza- 
tion over fiQ, the difference between x^ n an d the marginalized x 2 S n ^ s j us t a constant pC)| . 
We will adopt x 2 n as the goodness of fit between theoretical model and Snla data. 

The second set of data is the Baryon Acoustic Oscillations (BAO) data from SDSS 



DR7 [23 1, the datapoints we use are 



Hi. 2 



and 



do. 



r s (zd) 
D v (0.2) 

r s {zd) 



ZV(0.35)' 

where r s (z<i) is the comoving sound horizon at the baryon drag epoch and 



D v {z) 



dx 



H(x) J H(z) 



1/3 



encodes the visual distortion of a spherical object due to the non Euclidianity of a FRW 
spacetime. The inverse covariance matrix of BAO is 



^M,bao 



30124 -17227 \ 
-17227 86977 J 



The x 2 °f the BAO data is constructed as: 



Xbao 



Tsy-1 



Y 

M,bao ' 



where 



Y 



d .2 - 0.1905 
do.35 - 0.1097 



The CMB datapoints we will use are (R,l a ,z*) from WMAP7 [24]. z* is the redshift 



of recombination [32], R is the scaled distance to recombination 

dz 



R 



E z) 



and l a is the angular scale of the sound horizon at recombination 



la = K 



rja*) 



where r(z) = L dx j 'H(x) is the comoving distance and r s (a*) is the comoving sound 
horizon at recombination 

p c s (a) 
Jo a 2 H(a) 



-da, 



where the sound speed c s (a) = 1/y 3(1 + R\,a) and R^ = 3Q b /A£l\ is the photon-baryon 
energy density ratio. 
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parameter 






^bO 


best-fit+£;+£ 


„ 0t , 7 -U.UI6, -U.U26 
u - zo ' +0.009, +0.020 


U - DD/ +0.011, +0.021 


n n^/|- u - UU2 ' -U-0U4 

u - UO4 +0.001, +0.002 



Table 1: The best-fit values with la and 2a errors for f2 m o ; h and fibo i n the ghost dark energy 
model. 



The x 2 of the CMB data is constructed as: 



Xcmb — ^ ^M,cmb^> 



where 



X 



f l a - 302.09 \ 
R- 1.725 
z„ - 1091.3 



and the inverse covariance matrix 



^M,cmb 



( 2.305 29.698 -1.333 

29.698 6825.270 -113.180 
^ -1.333 -113.180 3.414 



The fourth set of observational data is 12 Hubble evolution data from [25] and p6| . 
Its xjj is defined as 



2 _ A [H(zj) - H ob (z t )} 2 

xh- 2_> ~2 • 

1=1 * 



Note that the redshift of these data falls in the region z € (0, 1.75). 

The last set we will use is the Big Bang Nucleosynthesis (BBN) data from 
whose x 2 is 

2 (n b0 h 2 - 0.022) 2 



Xbbn 



0.002 2 



In summary, we have 



/(.tot Xsn Xcmb Xbao Xh Xbbn 

and we assume uniform priors on all the parameters. 
3.3 Fitting Results 

The best-fit values and errors of parameters are summarized in Table |]. We also list the 
best-fit values of the corresponding parameters in Table ^ for comparison. The best-fit 
values of £l m o and h are slightly smaller than corresponding ones in the ACDM model. In 
Figure [3], we plot the ID marginalized distribution probability of each parameter. The 2D 
contour is plotted in Figure || from which we can see that there exists a strong correlation 
between f2 m o and h as we expected in subsection 3.1. 

With the best-fit value of O m o = 0.257, the transition between the dust phase and 

i 

the de Sitter phase occurs at a* = Ugs*-) 3 = 1.23. The universe begins to accelerate 
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parameter 




h 




best-fit 


0.273 


0.703 


0.045 



Table 2: The best-fit values for the ACDM model, using the same data sets. 
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0.048 0.050 0.052 0.054 0.056 0.058 



0.22 0.24 0.26 0.28 0.30 



Figure 4: ID marginalized distribution probability of h, £lbo and ft m o 
0.29 P ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' "~ 



0.28 



0.27 



0.26- 



0.25 - 



0.24- 



0.23 



0.22 



SN+Hubble+BBN+BAO+CMB 




0.62 



0.64 



0.66 



0.68 



0.70 



Figure 5: 68% and 95% contour plot in fi m o — h plane. The red dot in the figure stands for the 
best-fit value. 



at a acc — ~r — 0.615, or in terms of redshift, z a cc — 0.625. And the present EoS of DE 
u; = - TT -±— r = -0.796. 



f2m0+l 



X 2 of best-fit value of this model is Xmin = 607.192 for dof = 575. The reduced x 2 
equals to 1.056 which is acceptable. But Xmin ^ s l ar g er than the one for the ACDM model, 
Xacdm = 554.264. A similar conclusion is also reached by other authors using different 
data set [33]. 
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Figure 6: c 2 s ~ a where a* 



It is not hard to understand why in this model, Xmin 1S l ar g e > compared to the ACDM 
model. Recently many model independent studies on dynamics of DE show that current 



observational data favor cj ~ —1, at least at low redshift |34j]. But from (2.7) we can 
see that uq ~ — 1 requires a small O m o (^mo ~ 0), which goes against CMB and BAO 
observation. As a result, when we combine these different observational data sets to do 
joint likelihood analysis, the final \ 2 becomes large. 

4. Adiabatic Sound Speed and Interaction 

In this section, we will return to the simplified model introduced in section 2, and further 
study this model, where there are only two components, DE and CDM in a flat universe. 
The squared adiabatic sound speed of the ghost DE model is found to be 

2 P 11 

ct = -A- = , — < o. 

PDE 2 ( a_) +1 



The evolution behavior of the squared adiabatic sound speed is shown in Figure |6[ One 
can see from the figure that c 2 leaps from — | to at the time a*. Before a*, c 2 is less than 
zero, but after a* ~ 1, c 2 is approximate to zero. However, it is always negative! 

Negative squared adiabatic sound speed may cause some problems. For example, if 
a fluid evolves adiabatically and has no interaction with gravity or other fluids, negative 
squared sound speed implies an instability under perturbations. To see this clearly, let 
us consider a Newtonian argument. Making use of Euler equation, continuity equation of 
fluid and Poisson equation in Newtonian gravity, we have 



5 + 2H5 + c 2 s k 2 5 - 4irGp5 = S Q + S in , 

where 5 is energy density perturbation of the fluid. The terms on the r.h.s are source terms. 
The first term Sq represents interaction with other fluid, while the last term Si n stands for 
intrinsic entropy perturbation in the fluid itself. The 4irG term on the l.h.s represents the 
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effect of gravity. If there are no source terms, the negative c 2 term leads to an instability. 
But the presence of source terms will change 8's behavior and makes the thing complicated. 

In the ghost DE model, DE interacts with gravity only and there is no Sq or Si n term. 
Thus the squared sound speed criterion mentioned above shows DE will be unstable under 
perturbation in Newtonian limit, where we have neglected relativistic effects. However in 
a full relativistic treatment to discuss the stability of this DE model under linear perturba- 
tion, we need to define gauge invariant variables and solve perturbed Einstein's equation 
and conservation equations. However, it is beyond the scope of this note. 

From the above discussion, we see that negative squared adiabatic sound speed may 
cause a potential instability. But it can be improved if the source terms are present j39| . 
Therefore in the rest of this section, we introduce the direct interaction between DE and 
CDM and study the evolution dynamics of the model. 

In this case, Friedman equation still reads 

H 2 = — (aH + p m ) 

and conservation equations are modified to be 

Pm + 3Hp m = Q, 
Pde + 3(1 + u)Hp DE = -Q, 

where Q denotes the interaction between DE and CDM. Since a generic form of Q is 
not available, we consider three forms which are often discussed in the literature: Q = 
3aHpDE,3/3Hp m and 3jHp tot , where a, ft, and 7 are three constants, (see e.g. for 
more references). These forms imply that energy transfers in a Hubble time is proportional 
to energy density of DE, CDM and DE+CDM, respectively, and that energy transfers from 
DE to CDM if a, ft, 7 > 0, and vice versa. 

In terms of dimensionless quantities, we have 

i = n DE + n m , (4.1) 

= tl m + 2^0 m + 3HSl m , (4.2) 
-^Q = n DE + 2^n DE + 3(1 + oo)HQ DE , (4.3) 

where Qde = jfiiPDE, = jffiPm is dimensionless energy density of DE and CDM, 
respectively. In addition, by use of the linear relation between pjy E and H, we have 

Hn DE = H tt DE0 = const. (4.4) 

Combining fl4,l|) and (|4.2|) , one can get 

-fiflE + 2^ (1 - n DE ) + 3H(1- n DE ) = ^Q, (4.5) 
Putting ( |Q| ) to eliminate Q, we arrive at 

2H + 3loH 2 Q de + 3H 2 = 0. (4.6) 
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Then from ( [4.5[) and (4.4), one can have the equation of motion of Ode as 
2-n DE nTT ^ l-n DE 8?rG 



-n 



DE- 



DE 



where 



3a, when Q = 3aHpD E 



(4.7) 



n 



Q 



W when Q = 3 P r p™ 

3 ^II7^' when Q = SjHptot. 



Expressing this equation in terms of efolding-number N = In a, and making use of — -S^ 



^r-H tt DE0 /£l DE , we obtain 



, 2 - VLde „1- ^de 

""DE — — r o- 



n 2 



n 



n 



DE 



Q- 



Using flQ^flOp and (pj~7D , we can get the EoS of the DE 

1 2 n G 



2 — £Ide 3 2 — 



while the EoS of the total fluid is 



Utot 



DE 



1 2 H O 

-1 - = LdiLnf?, 

3ij2 



(4. 



(4.9) 



(4.10) 



in the second equality, we have used (|4.6| ). In addition, the deceleration parameter is given 
by 



da 1 — 20 



DE 



QMPE 

a 2 2-fl DE 2 -Q.de' 



and the squared speed of sound reads 

Pde 



q 



DE 



i 



DE 



2 — 0_d_B 



(4.11) 



(4.12) 



We will solve Qde, w, ujtot and analytically for each Q in the following. However, for 
the sake of briefness, we will discuss the case of Q = 3aHpDE only in some detail. 

4.1 Q = 3aHpD E 

In this case, (lOT) becomes 



Its solution is 



-q' de (2 - q de ) + 30^ - 3 (a + 1) q 2 de = o. 



3iV + C = 2lnn DE - 1 + 2a In 11 - (a + l)n DE \ 

1 + a 



(4.13) 



where the integration constant C = 21uQdeo — -jbf l n 1 1 — (" + 1) O^eoI; Od_eo is the 
dimensionless energy density of DE at present. 
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Figure 7: Qde ~ where N is efolding-number. Here £Ideo — 0.75. The blue, green and red 
curves correspond to a = 0, 0.1 and 0.2, respectively. 



If a < 0, from the solution ( 4.13| ), we can see that CIde can be larger than 1 at late time 



which is unphysical. This unphysical result comes from the fact that the assumption on 
energy transfer rate Q oc pde is oversimplified. For a < 0, DE will gain energy from CDM 
and energy density of CDM will become less and less. At some time, p m becomes zero, 
and it is impossible to keep on transferring energy to DE. But due to the oversimplified 
assumption on Q, CDM will continue to lose energy which is incorrect physically. Therefore, 
we will presume a > in this subsection. 

From ( |4. 13 ) one has (a + \)Q.de < !• It means that for a > 0, Qde will tend to 



a constant 1/(1 + a), rather than 1. The relation of £Ide ~ N is shown in Figure [7|, 
it also indicates that when a is larger, the evolution of £Ide will be flatter since more 
energy are injected into CDM. Of course there is a upper limit for a, as £Ide must be 
able to reach its present value Qdeo ~ 0.75. For a > 0, the coincidence problem can be 
alleviated excellently, and if a ~ — 1 + 1/Q,deo, this problem is completely solved. The 
similar situation occurs as Q ~ ptot (see section 4.3). 
In this case, the equation of state of DE is 

1 + 2q 

oj = — 



2 — O.DE 

We plot the relation uj ~ N in Figure ||. We have shown in section 2, it is impossible for 
uj to cross phantom divide without interaction. Nevertheless from this figure we can see 
that the situation is changed with the help of interaction term, oj will cross —1 from the 
quintessence regime to phantom regime and approach to —1 — 2a at late time when 
The equation of state of the total fluid is 

n „_x Qde 
utot = ~ (1 + 2a) 







DE 



and it is plotted in Figure [| As we expect again, uitot can be smaller than —1 in the presence 
of interaction term and reaches its asymptotic value —1 — 2a < —1 at late time. Therefore, 
in this model the universe will end with the big rip singularity in the future [36, 37 1. And 



-13- 




-1.2 



N 

Figure 8: w ~ N where N is efolding-number. Here £Ideo — 0.75. The blue, green and red curves 
correspond to a = 0, 0.1 and 0.2, respectively. 




N 

Figure 9: ujtot ~ -/V where iV is efolding-number. Here ftoEQ = 0.75. The blue, green and red 
curves correspond to a — 0, 0.1 and 0.2, respectively. 

from this figure, we can see that the universe begins to accelerate earlier with larger a. 
Note that according to the calculations in section 2 and section 3, the universe begins to 
accelerate at z acc ^=o = 0.6 in the case without interaction. Thus the acceleration of the 
universe occurs at z acc> a > 0.6 when the interaction is present if one keeps f2 m o = 0.257. 
Finally, the squared speed of sound reads 



P 



Pde 



-2(1 + 26) 



1 — Qde 



n 



DE 



\2 ' 



(4.14) 



the squared sound speed is always smaller than because we assumed a > 0. This means 
that this form of interaction cannot make (? s positive. 

4.2 Q = 3PHp m 

In this case, ( |Q| ) becomes 



-(2 - n DE )n' DE + 3(1 - p)n DE (i - n DE ) = o. 



-14- 




its analytical solution reads 

3(1- P)N + C = 21nO DB -ln(l-0 0£ ), 

where the integration constant C = 2\yiQdeq — hi(l — Qdeo)- We can see from the 
analytical solution that Qde varies from at early time to 1 at late time. The EoS of DE 
is 

1 2/ g 1 l-^DE 

2 — Qde 2 — Qde Qde 
One can see in this case that at early time w —> —(3 /Qde as Qde — > 0. The plot is shown 



in Figure 10. If f3 > 0, ui will increase from — oo to some local maximum at some point 
and then decrease to its asymptotic value —1 at late time. Unlike the situation we have 
discussed in subsection 4.1, uj will cross —1 from phantom regime to quintessence regime 
in this case. If fi < 0, u will decay monotonically from +oo to —1 and never cross — 1. 
However, no matter what the value of (3 is, the late-time asymptotic behaviors of u are all 
the same. 

The EoS of the total fluid is 

2/3 + (l - 2/3) Qde _ \ -/3 at early time 



Utot 



Qde I —1 at late time. 



The reasonable value of /3 should be |^| < 1. Thus the big rip singularity will be avoided 
in this case. 

Finally, we give the squared speed of sound 



1 - Q 



DE 



(2-Q 



DE 



4 - 5Q DE + 2Q 2 DE 



(2-Q 



DE 



Q 



DE 



from which we can find that 



-> —2(3 /Qde at early time. Therefore, this interaction 
term will not make (? s positive as well if /3 > 0. Furthermore, at early time the speed of 
sound will be larger than speed of light if ft < 0. 
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4.3 Q = 3^Hp tot 

In this case ([Q|) becomes 



. 2 - Q DE 1 - &DE 
' n DE—Fv2 + 3 r. = 3 ^ 



fiflg fiflE 

Its analytical solution reads 

3N + C = -^lnn DE - (1^1) ln|l-7-0 DB |, 



1 — 7 \ 1 — 7 

where C = ^ In f} D £ - In |1 - 7 - Qdeq\ ■ The EoS of DE 

1 7 1 

w = -- ^ 2- 



2 — £Ide 2 — Qde 
Once again w will diverge at early time when the interaction is present. The EoS of the 
total fluid is 

27 Ude 

Utot 



2 — &DE 2 — Q.DE 

and the squared speed of sound reads 

2 o 1 — ^de . - 3fi£)£; — 4 

C s = ~ 2 2 + 2 ^ 2 

(2 — O^s) (2 — f2_DE) 

which is also divergent at early time £Ide 0. 
5. Conclusion and Discussion 

In this note we investigated a DE model whose energy density is proportional to Hubble 
parameter with a coefficient which is roughly order of Aq CD . It gives the right order of 
magnitude of observed energy density of DE. We studied its cosmological evolution. In 
this DE model, the universe has a de Sitter phase at late time and begins to accelerate at 
redshift around z acc ~ 0.6. 

We also fitted this model with observational data including Snla, BAO, CMB, BBN 
and Hubble parameter data. The best-fit values of parameters of the model are f2 m o = 
0.257, h = 0.662, Q^q = 0.054. However, the minimal x 2 gives Xmin = 607.192, while in the 
ACDM model, Xacdm = 554.264 for the same data sets. Namely the simple x 2 analysis 
seemingly implies that current data do not favor the ghost DE model, compared to the 
ACDM model. Clearly this result is not conclusive, further study is needed. 

We also found that the squared sound speed of the DE is negative, which may give 
rise to a potential instability of the model under perturbation. We further studied the 
cosmological dynamics of the model by considering there exists some interaction between 
DE and CDM. Three kinds of interaction forms are discussed. In all cases, the negative 
squared sound speed is still there with a proper coefficient for the interaction terms. Clearly 
the potential instability should be studied seriously by investigating linearized Einstein 
equations, not just calculating the squared sound speed of the DE, which is currently 
under investigation. If the instability indeed exists, then the ghost DE model has to be 
further modified or to be abandoned. 
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